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We present a complete analysis ol threshold resummation effects on direct light and heavy gaugino 
pair production at the Tevatron and the LHC. Based on a new perturbative calculation at next- 
to-leading order of SUSY-QCD, which includes also squark mixing effects, we resum soft gluon 
. . . radiation in the threshold region at leading and next-to-leading logarithmic accuracy, retaining at 

fy^ , the same time the full SUSY-QCD corrections in the finite coefficient function. This aUows us to 

^~~^ • correctly match the resummed to the perturbative cross section. Universal subleading logarithms 

^^ ' are resummed in full matrix form. We find that threshold resummation slightly increases and 

Cn , considerably stabilizes the invariant mass spectra and total cross sections with respect to the next- 

^.,^. to-leading order calculation. For future reference, we present total cross sections and their theoretical 

j^ ' errors in tabular form for several commonly used SUSY benchmark points, gaugino pairs, and hadron 

^rj ^ collider energies. 

^— I . 

I \'. I. INTRODUCTION 

^■ 

I The Minimal Supersymmetric Standard Model (MSSM) [1] continues to be one of the best-motivated and most 

Qh, intensely studied extensions of the Standard Model (SM) of particle physics. It introduces a symmetry between 

^ ■ fermionic and bosonic degrees of freedom in nature and predicts one fermionic (bosonic) supersymmetric (SUSY) 

partner for each bosonic (fermionic) SM particle. Consequently, it allows to stabilize the gap between the electroweak 

and the Planck scale [2] and to unify the three gauge couplings at energies of 0(10^^) GeV Q. It also contains a stable 

lightest supersymmetric particle (LSP), which interacts only weakly and represents therefore an excellent candidate 

^ , for cold dark matter [Jl. Spin partners of the SM particles have not yet been observed, and in order to remain a 

viable solution to the hierarchy problem, SUSY must be broken at low energy via soft terms in the Lagrangian. As a 

^v ' consequence, the SUSY particles must be massive in comparison to their SM counterparts, and the Tevatron and the 

fvq LHC will perform a conclusive search covering a wide range of masses up to the TeV scale. 

The production of SUSY particles at hadron colliders has been studied at leading order (LO) of perturbative QCD 
since the 1980s 0. More recently, previously neglected electroweak contributions Q, polarization effects [7|, and 
the violation of flavor [8] and CP symmetry \9\ have been considered at this order. Next-to-leading order (NLO) 
corrections have been computed within QCD since the late 1990s [lOl \1)\ and recently also within the electroweak 
theory [12] ■ Resummation at the next-to-lea ding logarithmic (NLL) level has been achieved in the small transverse- 
momentum region for sleptons and gauginos [ISl, in the threshold region for sleptons, squarks and gluinos 14], and 
for sleptons and additional neutral gauge bosons also in both regions simultaneously [15] . 
5— i As the fermionic partners of the electroweak gauge and Higgs bosons, the four neutralino (x^, i = 1, ■ • ■ ,4) and 

. . 1 chargino (x^ , i = 1,2) mass eigenstates are superpositions of the neutral and charged gaugino and higgsino interaction 
eigenstates. Their decays into leptons and missing transverse energy, carried away by the Xi LSP, are easily identifiable 
at hadron colliders. The lighter mass eigenstates are accessible not only at the LHC with center-of-mass energies y/S 
of 7 to 14 TeV, but also at Run II of the Tevatron {y/S = 1.96 TeV), where the production of Xi X2 pairs decaying into 
trilepton final states is a gold-plated SUSY discovery channel [la] . For this particular channel, threshold resummation 
has been studied previously within the minimal supergravity (mSUGRA) model [17| . It was found that threshold 
resummation could increase the total cross section by up to 4.7% at the Tevatron relative to the NLO prediction, i.e. 
was significant even relatively far from the hadronic production threshold. 

In this work, we extend and improve on these published results in several respects. First, we extend the original 
NLO calculation [ll| by also including squark mixing effects. Second, we include not only the QCD, but also the 
SUSY-QCD virtual loop contributions in the hard coefficient function of the resummed cross section, which therefore 
reproduces, when expanded, the correct NLO SUSY-QCD cross section in the threshold region. Third, we resum 



^ 



o 



o 
o 



X 



IklasenQlpsc . in2p3T&] 




Xi 9 




Xi 1 




I 



Xi 



q ^ ^Xj q — ^^""'^^Xj q ^ ^^""'^^Xj 

FIG. 1: Tree-level Feynman diagrams for the production of gaugino pairs. 



not only the diagonal, but the full matrix contributions coming from the anomalous dimension, thereby including 
all universal subleading terms and full singlet mixing. However, we refrain from resumming constant terms that are 
known to factorize and exponentiate in Drell-Yan like processes, but not in more complex processes like gaugino pair 
production with several interfering Born diagrams. For the Tevatron, we consider not only the production of x-^ X2j 
but also of X2X2 ^^^ Xi Xi pairs. In particular the latter can have significantly larger cross sections than trilepton 
production due to the s-channel exchange of massless photons. For the LHC, we concentrate on predictions for its 
initial center-of-mass energy of v^ —1 TeV and include also the production of heavy gaugino (xi > X3,4) combinations, 
where threshold effects and direct gaugino pair production (as opposed to the production from squark and gluino 
cascade decays) will be more important. However, we will also show cross sections for \fS =14 TeV for comparison. 
The remainder of this paper is organized as follows: In Sec. HIl we present briefly our LO and NLO SUSY-QCD 
calculations, focusing on squark mixing, the ultraviolet renormalization procedure and the dipole subtraction method 
employed for the cancellation of infrared divergences among virtual and real contributions. In Sec. lIIIl we describe the 
threshold resummation formalism in Mellin space, starting from the refactorization of the perturbative cross section 
and the exponentiation of the eikonal function, giving then the explicit form of the resummed logarithms at NLL 
order, and finally describing various ways to improve on the original formalism. Sec. IIVI contains numerical results for 
various gaugino pair production cross sections at the Tevatron and at the LHC in both graphical and tabular form. It 
also includes a comprehensive analysis of theoretical errors coming from scale and parton density uncertainties. We 
summarize our results in Sec. [V] Our conventions for the couplings of quarks and squarks to weak gauge bosons and 
gauginos are defined in the Appendix. 



II. GAUGINO PAIR PRODUCTION AT FIXED ORDER IN PERTURBATIVE QCD 



At hadron colliders, two gauginos Xi ,j with masses ra^^ . and four-momenta pi,2 can be produced at LO of pertur- 
bation theory only through the annihilation of quarks q and antiquarks g' with four-momenta Pa,&, 



q{va) q'ipb) -^ Xi{pi) X]{P2), 



(1) 



since gluons couple neither to electroweak gauge/Higgs bosons nor to their SUSY partners [43] . As shown in Fig. [l] 
the production can proceed either through the s-channel exchange of a photon (7), neutral {Z) or charged (W) weak 
gauge boson, depending on the charge of the final state, or through the t- and w-channel exchange of a squark (q), 
where the Mandclstam variables 



(Pa+Pb)'^ , t = (Pa-Pl)^ , and U = {Pa-P2f 



and their mass-subtracted counterparts t ^^^ ^ — t — m| . . and u^^- ■ 



u — m^, . are defined in the usual way. 



(2) 



A. Leading order cross section 

The LO differential cross section for given (anti-)quark helicities ha.b can be expressed in the compact form 
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TABLE I: s-channel charges for gaugino pair production. The masses of the W- and Z-bosons are mw and mz- 
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by employing the generalized charges 
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where the s-channel charges Q^c ^''^d Q3f are given in Tab. HI rriq is the mass of the squark exchanged in the t- or 
M-channel, X,Y ^ {L,R} relate to the left- and right-handed (s)quark couplings to weak gauge bosons and gauginos 
listed in App. |Al and the couplings of the latter among each other depend on the neutral and charged gaugino- higgsino 
mixing matrices N, U and V through the bilinear combinations [l[ 
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a = e^/(47r) is the electromagnetic fine structure constant, xw ~ ^ — c^ — s^ — sin^ 6w is the squared sine of the 
electroweak mixing angle, and Ca — 3, Cf = 4/3 and as = g^/(47r) are the QCD color factors and coupling constant. 
Unpolarized cross sections are obtained by averaging the polarized ones with 
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B. Virtual corrections 

At NLO of SUSY-QCD, ©(cts), the cross section for gaugino pair production receives contributions from the 
interference of the virtual one- loop diagrams shown in Figs. [^HHwith the tree-level diagrams shown in Fig. [l]on the 
one hand and from real gluon (Fig. [5]) and (anti-)quark emission diagrams on the other hand, where the latter are 
obtained by crossing the final-state gluon in Fig. [5] with the initial-state antiquark (Fig. [6]) or quark (not shown). All 
diagrams have been evaluated analytically with self-written FORM programs and cross-checked independently with 
self-written MATHEMATICA programs. 

The virtual self-energy diagrams for left- and right-handed quarks qL^R = Pl.rQ = 1(1 T 75)9 (Fig- HI third line). 
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FIG. 2: Self-energy insertions (top) and contributions (bottom) to the production of gaugino pairs. 

expanded as usual into vector (V) and scalar (S) parts S(p) = [S)^(p^)|> + S£(p^)]Fl + {L ^ R), as well as those for 
squarks (Fig. [2l fourth line), 
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contain ultraviolet (UV) divergences in the scalar integrals i?o,i,...(Pj"^i,''7i2) [3) which exhibit themselves as 1/e 
poles in £) = 4 — 2e dimensions [4J] . They must therefore be absorbed through a suitable renormalization procedure 
into the fundamental wave functions, mass parameters, and coupling constants of the SUSY-QCD Lagrangian 
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(15) 



The two components of the unrenormalized squark field qP correspond originally to the left- and right-handed chiralities 
of the unrenormalized SM quark field g", but mix due to the fact that soft SUSY-breaking and Higgs terms render 
the 2 X 2-dimensional mass matrix (to?)° non-diagonal [l[. In Eq. ([T5|) we have diagonalized this mass matrix with 
the squark rotation matrix i?^°, so that the components i = 1 (2) of the squark field correspond to the squark 
mass eigenvalues m?. . The squark self-energies in Eqs. (fT2|) - ([T4|) thus also carry indices i,j — 1,2 corresponding to 
the (outgoing and incoming) squark mass eigenstates. Multiplicative renormalization is achieved perturbatively by 
expanding the renormalization constants. 
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FIG. 3: Vertex correction insertions (top) and contributions (bottom) to the production of gaugino pairs, 
with the usual factor of 1/2 for the (s)quark wave functions. The renormahzed self-energies are then 
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for squarks. 

We choose to renormalize the wave functions in the MS-scheme, so that the definition of the quark fields corresponds 
to the one employed in the parton densities in the external hadrons [45|. In this scheme, the quark wave function 
counterterm 



5Z„ = SZ^f^+SZ^f) ^ith 5Zlf) = 5Z'f^ = -%^A and A = --7£; + ln47r, 



167r2 



(20) 



defined as the UV-divergent plus universal finite parts of the on-shcU counterterm — I]^^(rn^) — m^[I]J^'^(m^) + 
S^'^(to^)] — TOq[E|'jj.(m^) + S;|'^(to^)] [1^, is hermitian {5Zq = 5Z'^) and the same for left- and right-handed quarks. 
The superscripts g and g label the gluon and gluino exchange contributions, respectively, and ^e is the Euler constant. 
The squark wave function counterterms 
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FIG. 4: Box diagrams contributing to the production of gaugino pairs at NLO 
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for i^ j, defined similarly as the UV-divergent plus universal finite parts of the on-shell counterterms — ReE^j(TO? ) 

for i = i and 2 Re Ey (?Ti? )/(to? — fnjj) for i ^ j |jL2|, enter only through the renormalization of the squark mixing 
matrix, 
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since in the MS-scheme the gluon and gluino contributions for i = j in Eq. ([22) cancel each other. In the last step 
of Eq. (l26l) , we have made use of the fact that in the MS-scheme the squark wave- function renormalization constants 
are anti-hermitian matrices {SZq^j = —SZ~,^). 

The (s)quark masses are renormalized in the on-shell scheme to make them correspond to the physical masses. The 
quark mass counterterm is then defined by mgS^(rnq) -|- S'^(mq) [13] with the result 
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For our numerical results, we will set the masses of external quarks to zero in accordance with the collinear factorization 
of quarks in hadrons. The squark mass counterterm is defined by Re[Eii(m|.)]. The result is 
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FIG. 5: Gluon emission diagrams contributing to the production of gaugino pairs at NLO. 



Supersymmetric Ward identities link the quark-quark-gauge boson and quark-squark-gaugino vertices to the weak 
gauge-boson and gaugino self-energies. As the latter do not receive strong corrections at NLO, the former require 
no further renormalization beyond the one for the (s)quark wave functions discussed above. However, the artificial 
breaking of supersymmetry by the mismatch of two gaugino and {D — 2) transverse vector degrees of freedom must 
be compensated by a finite counterterm g = g[l — asCp/{STT)], effectively shifting the quark-squark-gaugino scalar 
coupling constant g with respect to the weak gauge coupling constant g 
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C. Real corrections 
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Apart from the (now UV-finite) virtual corrections da^^^ ' to the LO cross section da]^f^' described above, the NLO 
cross section 
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also receives contributions da^^ ' from real gluon (Fig. [SJ , quark (Fig. [S]) and antiquark (not shown) emission dia- 
grams, where the emitted parton carries four- momentum p^. In the Catani-Seymour dipole formalism |2l| . the real 
contributions are rendered infrared (IR) finite by subtracting from them their soft and collinear limits {pa,b ■ Ps ^ 0) 



'i<^abiPa,Pb) '^dVdipolc = ^ ^V°-^'''{pi,P2,P3;Pa,Pb) Ff\pi,p2\Pa3,Pb) + (« ^ 



(35) 



before integration over the three-particle final-state phase space. They can then be evaluated in four dimensions (i.e. 
with e — 0). In the case at hand of two initial-state partons and no colored final state particles at LO, the only dipole 
contribution comes from an initial state emitter and an initial state spectator, e.g. 
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FIG. 6: Quark emission diagrams contributing to the production of gaugino pairs at NLO. 



(see Eq. (5.136) of Ref. ^). The color charges Ta3,b and sphtting functions V^^^ (Eqs. (5J.45)-(5.148) of Rcf. 
[211 ) act on Born-hkc squared matrix elements, which are written here in terms of vectors |1, 2; a3, 6)2,06 in color and 
helicity space. These matrix elements involve an initial-state parton a3 with momentum parallel to Pa, 



P^3 = ^3,abPa , whcrC Xs^ab 
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and rescaled four-momenta of the final-state gauginos 
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where K^ = Pa +Pfc "^3 ^'^'i ^^ — PaS +Pfc ■ The phase space function Fj {pi,p2]Pa3,Pb) tends to zero with pa • ps 
and ensures therefore that the LO cross section is IR-finite. To compensate for the subtracted auxiliary dipole term 
'^'^ab {Po-TPb) ^ dVdipoioj the latter must be integrated analytically over the full phase space of the emitted parton, 
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and added to the virtual cross section. The integrated dipole term is defined explicitly in Eq. (10.15) of Ref. [2i|; it 
contains all the simple and double poles in e necessary to cancel the IR singularities in dcr^j^ . The insertion operators 
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FIG. 7: Associated production of a gaugino and a virtual squark, decaying subsequently into a gaugino and a quark. 



where /3o = 11Ca/6 - 2NfTR/3 and /3i = (17C^ - SC^iV/ - 3CFNf)/6 are the one- and two-loop coefficients of 
the QCD beta-function, Cp — 4/3, Tfi — 1/2, Ca — 3, and Nj is the number of quark flavors. They cancel the 
dependence of the hadronic cross section on the factorization scale ^f up to NLO accuracy. The insertion operators 
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depend on the factorization scheme through the term K^'^ (x), which vanishes in the MS-scheme, and also on the 
regular parts of the Altarelli-Parisi splitting distributions given by P^^ (x) = P°'^{x), if a ^ 6, and otherwise by 



i^gg(Xj 



Cf (l + x) and P39{x) = 2Ca 



1-x 
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(52) 



The last line in Eq. ([M)) contains therefore the finite remainders that are left after the factorization of collinear 
initial-state singularities into the parton densities in the MS-scheme at the factorization scale ^f- As guaranteed by 
the Kinoshita-Lee-Nauenberg and factorization theorems, the total NLO cross section is then not only UV-, but also 
IR-finite. 

Finally, one subtlety must still be addressed: in Fig. ^ the center and right diagrams of lines one and three proceed 
through a squark propagator, which can become on-shell if rriq > m^ and s > {rriq + m^Y' . To avoid double counting, 
the resonance contribution 



dcr'f == da{qg -^ xq) BR(g ^ xq) 



(53) 



must be subtracted from the gaugino pair production process using the narrow-width approximation, as it is identified 
experimentally as the associated production of a gaugino and a squark (Fig. [T]), followed by the decay of the squark 
into a gaugino and a quark (Fig. |8]). 
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FIG. 8: Tree-level diagram for a squark decaying into a gaugino and a quark. 

III. THRESHOLD RESUMMATION AT NEXT-TO-LEADING LOGARITHMIC ACCURACY 

In the previous section, we have demonstrated that all soft and coUinear (IR) singularities in the partonic NLO 
cross section dcr^j, either cancel among virtual and real corrections or can be absorbed at the factorization scale fip 
into parton density functions (PDFs) fa,b/A,B{xa,bTy''F), which represent probability distributions for initial partons 
a, h with longitudinal momentum fraction Xa^ in the external hadrons A, B (e.g. protons or anti-protons). The QCD 
factorization theorem guarantees that the observable hadronic cross section 

M^ (r) = ^ / dXadXbdz[Xafa/A{Xa,^J''^)][xbfb/B{xb,^^'^)][z<i(Tab{z,M'^,fI'^)]S{T - XaXbz) (54) 

ab 

can be obtained by convolving the process-dependent partonic cross section 

oo 

daabiz^M^f,-") = Y,a:{f,^)d<jil\z,M\fi') (55) 

with universal PDFs. For the sake of simplicity, we will identify in the remainder of this paper the factorization scale 
fip and the renormalization scale hr, which enters the partonic cross section through the perturbative expansion in the 
(reduced) strong coupling constant as(/i|j) — as{fi%)/{2Tr) and also explicitly beyond LO, with the common scale ji. 
M"^ represents the invariant mass squared of the (colorless) gaugino pair produced at LO, which carries the fractions 
T = M"^ / S and z = M"^ / s of the hadronic and partonic center-of-mass energies S and s — XaXbS, respectively. At 
LO, where no additional partons are produced besides the gaugino pair, the partonic cross section 

dai°)(z,M2,M2) = a^^XM^)5{l - z) (56) 

is independent of /i and has its support entirely at the point z — 1. At higher order in QCD, the cancellation of soft 
and coUinear parton emission among virtual and real corrections is restricted by the phase space boundary of the 
latter. This leads to logarithmic contributions a"[ln™(l — z)/(l — z)]_|_ with m <2n—l (see Eqs. (001) and (05])), which 
become large close to the partonic threshold at z — > 1 . Consequently, they spoil the convergence of the perturbative 
series and must be resummed to all orders in a^. 
By applying a Mellin transform 

F{N) = f dyy^-'F{y) (57) 

to the quantities F — <7ab, crab, and fa,b/A.B with y — t, z, and Xa.b in Eq. (j54|) . the hadronic cross section can be 
written as a simple product 

M'^^{N-1) = J2f-MN,^^')h/B{N,^l^)aab{N,M^fl'), (58) 

ab 

the large logarithms in z — )■ 1 turn into large logarithms of the Mellin variable A^, 

'^'^"'(^-')\ -ln™+iiV etc., (59) 



l-z ,^ 
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and it is thus possible, in the N —i' oo hmit, to retain only the leading powers in N. In particular, one can neglect 
among the TV- moments of the Altarelli-Parisi splitting functions at 0{as) 



P''{N) = 1^'Xn) - C, 



P'^iN) = l!,l\N) = Cf 



2 N{N 

2 + N + N^ 



1 ^ 1 

r + 1) ^k 



P'^iN) = jil\N) 



Tn 



N{N^ - 1) 
2 + N + N^ 



fc=i 

N ' 



Cf{ --2\nN 



^1^1' 



P''iN) = J^'KN) = 2Ca 



N{N + 1){N + 2) 
1 



Tr 

N 



and 



N{N-1) {N + l){N + 2) ^k 



^ 1 
2^ fc 



/3o 



(60) 

(61) 
(62) 



2CAlniV + /3o + 0(-) (63) 



those which are not diagonal. Here, we have introduced the variable N — iVe'*'^, where 7^ is the Euler constant. 

A. Refactorization 

The resummation of the large logarithms in Eqs. (|40)) and (|45|) is based on the observation that the separation of 
the non-perturbative PDFs fa,b/A.B ^^d the perturbative partonic cross section Uab in Eq. (|58p is not unambiguously 
defined. It is in particular possible to re-factorize the partonic cross section [22] 



2 ,,2\ 



aab{N,M\^i 



i^.(M-..-) V^^^!:;t-er^ ^..(^,M^)+^(^ 



fa,a{N,^l')h,b{N,^i^ 



(64) 



into a hard function Hah, which is non-singular in and in fact independent of N ^ a ratio of PDFs in partons ^a.b and 
fa,b defined at fixed energy and longitudinal momentum fraction, respectively, and a function 5a6, which describes 
the emission of soft giuons and can thus be computed in the eikonal approximation. The hard function 



H^b{M\^,') ^ ^a:(^^)i?(:)(Af^/.2) 



n=0 



can be expanded as a power series in a^. Its LO and NLO coefficients read 

H^XM^^?) = aS)(M2)and 



r(i)^A/r2 ,,2 



(0), 



H^^>{M\^,^) = aZ'{M^) 



A. 



o+(7«+7^^0 



In- 



M 



#2 



where Aq represents the IR-finite part of the renormalizcd virtual correction 



aiP(M^;.2) = a. 



Atth^Y r(l-e) fA-2 A-i 



M2 J r(l - 2e) V e 



(0), 



+ l-i+A <Tl",^(Af^) + 0(e). 



The PDFs satisfy the evolution equations [23 



dfc/aiN,n'^) -^ ^ 2^^ , f .r 2^ 
Q^^2 = 2^1bc{N,asifi ))fb/a{N,^l ), 



(65) 

(66) 

(67) 

(68) 
(69) 



governed by the (gauge-independent) anomalous dimensions jaa'^N^aglp"^)) = J2n^^sit^'^)laa'(^) of the composite 
Wilson operators with 7^ J (N) given above, and [23] 



91n/i^ 



7a(as(M ))^a/a{N,^J- ), 



(70) 



where the (gauge-dependent) anomalous dimensions ja — ^/Za dZa/dlnfi'^ — J2n o"(M^)7a of the fields a correspond 
in the axial gauge Q to the TV- independent (virtual) parts of "faaiN) [4^. In the singlet /non-singlet basis and to 
one-loop order, the evolution equations for f^/^ can be solved and written in the closed exponential form [25j 



f,/a{N,fl^) =. J2^bc{N,^l^^ll)h/a{N,^ll), 



(71) 
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where the evohition operator E{N , ji^ , jif^) satisfies the same evolution equation, Eq. ((69|) . as fc/a{N, fi^). The ratio 



fa/a[N,fJ.^) 



(72) 



has been shown [22[ to exponentiate to all orders into a gauge-independent, but scheme- and scale-dependent function 



lnA,(iV,M2,M') = / d^^ / 

Jo 1-z J(i 



dg^ 



'(1-^)'"=^'^ 9 



^„(a,(g^))-i?„(a,((l-zr=M^)) 



(73) 



with TTis = and consequently Ba = in the MS-scheme, which describes the soft and collinear gluon radiation from 
the initial partons, and a gauge-dependent, but scheme- and scale independent function 



\nUa{N,M^) = 
which can be combined with the soft function 



dz 

1-2 



dz 

1-2 



1 



-iya{as{{l-z)^M^)), 



-Xa{a,{{l - zfM^))Sab 



into the gauge-, scheme- and scale-independent function 



\nAabiN,M^ 



dz 

1-2 



-Dabias{{l-z)Hl^)). 



(74) 



(75) 



(76) 



The functions Ua and Sab can be computed in the eikonal approximation. They depend on the cusp anomalous 
dimension i^a = 2Ca(7Coth7 — 1) [26| of the considered process and describe soft wide-angle radiation. 



B. NLL approximation 



The coefficients 



b 



Aaias) = £a;M("), Ba{a,) = £a^Bi"\ and D,,(a,) = £ a^i^i' 

n— 1 n— 1 

can be expanded as power series in a^, and the NLL results read (221. |27J 

[Vl8 6 J 9 ^ 

B^i) = (in the MS scheme), and 
dW = (since ^W=AW) 



(77) 

(78) 
(79) 

(80) 
(81) 



with vi — AC a and Ca = Cp and Ca for an incoming quark (a = q) and gluon {a — g), respectively. Note, however, 
that D^j' ^ [2^. After the integrations in Eqs. ([75)1 and ([75)1 have been performed, the partonic cross section in 
Eq. (|64p can be written in the closed exponential form [28| 



aab{N,M^,fi^) = Hab{M^,f^^)exp[gab{N,M^,fi^)]+0[- 



N 



(82) 



Here, the perturbative coefficients of the hard function 



nab{M',f^') = J2a^{^^}nl:HM',^,') 



(83) 



n=0 
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have been redefined with respect to those in Eq. (|65|) in order to absorb the non-logarithmic terms resulting from the 
integration, i.e. 



(85) 



H^:,'{^I^,^.'') = a:,^(Af^)and (84) 

The coefficient function T-L)^^ {M'^ , n"^) given here agrees with the one presented in Eq. (115) of Ref. [17| except for 
their last three terms. While their last term corresponds to the flavor-diagonal collinear improvement to be discussed 
below, the two other terms represent leading and next-to-leading logarithms and should therefore not be present. 
Furthermore, the two terms in (ln47r/x^/(5^ — 7b) should be squared individually, not together, and the virtual 
correction My defined in Eq. (116) of Ref. [I4I should include the complete SUSY-QCD contributions and not 
only their UV-singular parts. The function Qab takes the form 

gab{N,M\^?) = Lgi'^{\)+gl^,\\,My^i^) + aJ^l\x,Myfi^) + . . . (86) 

with X = asPoL and L ^\nN. The first term in Eq. (1551) . 

LgilH^) = -^(A«+4'^)[2A+(l-2A)ln(l-2A)], (87) 

collects the leading logarithmic (LL) large- iV contributions L{asL)'^ and depends on A\ I only. The coefficients A^ I, 
A^^ I and D^j, determine the function 

2/3o.9f,)(A, AfV^') = {Ai'^ + 4'^) ln(l - 2A) In ^ 

+ (AW +ylW)§[2A + ln(l - 2A) + ^^^(1 - 2A)] 

- (4^) + 4'^)^ [2A + ffi(l - 2A)] + i?« ln(l - 2A), 
Po 

which resums the next-to-leading logarithmic (NLL) terms (a^L)". Similarly, the functions gJl resum the N"LL 
terms and depend on the coefficients A^^j^ ' , A^^'l and D^^^ with 1 < k < n. 

C. Collinear improvement 

Up to this point, we have systematically neglected all terms of 0{1/N). However, since the dominant 1/A^-terms, i.e. 
those of the form a"L'^'^~^/N, stem from the universal collinear radiation of initial state partons, they are expected to 
exponentiate as well. This has been proven to next-to-next-to-leading order for deep-inelastic scattering and Drell-Yan 
type processes [23] and can be achieved by making the replacement (cf. Eq. ([U7)) below) 

n[\'^Hl^,'+L^^-^U^:,\ (89) 

i.e. by including the corresponding subleading terms of the diagonal splitting functions jaa,bb{N) in Eqs. (|60p and 
(1631) . Carrying on with this argument, it is even possible to resum the terms of 0{1/N) coming from the diagonal 
and non-diagonal splitting functions by identifying the terms [30| 

^5ll'(A) = -^(AW+AW)[(-2A)ln(l-2A)]+... (90) 

with the LL approximation of the QCD evolution operators Eat defined in Eq. (j7ip and then promoting the LL to 
the full one-loop approximation E^jJ . The resummed cross section, Eq. (|82l) . can then be written in a coUinearly 
improved form as 

aabiN,M',^i') = Y,4'J'iN,MyN',^l%)na'b'{M^,^^l) exp[ga'b'{N,M^til)]Eil}{N,MyN^,fil), (91) 

a'b' 
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where the dependences on the factorization and renormahzation scales have been recovered exphcitly, the colhnearly 
improved hard coefhcient function TLab is expanded as usual as a power series in Usip'^) and its LO and NLO coefficients 
read now 

n^:^{M^^^') = ai°)(M2)and (92) 



?^«(M^^^) = a^iM^) 



A + y(AW+A(^ 



b 



(93) 



and the Sudakov exponential function Qab is expanded in the same way as Gab in Ec^. (I86p with 



L~9il\>^) = :^(AW+4^))[2A + ln(l-2A)] and (94) 



2/3o5i'^(A,MVA^') = (AW+4'))[2A + ln(l-2A)]ln^ 

+ (A(i) + 4'^)^ [2A + ln(l - 2A) + i ln2(l - 2A)] 

- (A(2)+4'))l[2A + ln(l-2A)] 
Po 

+ (-27W-27l') + i?W)ln(l-2A). (95) 

Here, the anomalous dimensions 7^ I in g^i, {X, IVP /ji^) have been introduced to cancel the NLL terms in the one- loop 

approximation E^f^ of the evolution operators. 

For the Drell-Yan process, it has been suggested that also the constant terms in the hard coefficient function 
"^ab (-^^7 /^^) can be exponentiated, since these terms factorize the complete Born cross section, include finite remain- 
ders of the infrared singularities in the virtual corrections and are thus related to the corresponding singularities in 
the real corrections giving rise to the large logarithms [3l|. While this choice is supported by an explicit two-loop 
calculation [S^] and can be applied to other Drell-Yan like processes [Ij], gaugino pair production does not proceed 
through a single s-channel diagram, and the virtual corrections thus factorize only at the level of amplitudes, but not 
the full cross section. Resumming some or all of the finite terms in the hard coefficient function Ti^^ {M"^, /x^) into an 
exponential as in Ref. [iJl seems therefore not to be justified in this case. 

D. Matching and inverse Mellin transform 

As mentioned above, the large logarithms, which spoil the convergence of the perturbative series and must be 
resummed to all orders, appear close to production threshold. Conversely, the perturbative cross section should be 
valid far from this threshold. To obtain a reliable prediction in all kinematic regions, both results must be consistently 
matched through 

(res.) (f-o.) (cxp.) ^nc\ 

Tab = <b '+<^ab -<b ' (96) 

i.e. by subtracting from the sum of the resummed (res.) cross section in Eq. (j9ip and the fixed order (f.o.) cross 
section in Eq. (|34)) their overlap. The latter can be obtained by expanding (exp.) the resummed cross section to the 
same fixed order as the perturbative result. At ©(a^), one then obtains 

ai^■)(iV,M^^^) = Hi°)(M^^^)+a.H«(A/^^^) 

2i - In ^) 5: bi\\N)n'^\M^ ^2) + 7ii°)(^f^ ^.'hil\N)] 

- aM^M\ ^?) [L\A^) + 4^)) - 2L(7i^) + 7^^^)] ■ (97) 

After the resummed result and its perturbative expansion have been obtained in Mellin A^-space and multiplied 
with the A^- moments of the PDFs according to Eq. (|58p , an inverse Mellin transform 



"^St^w - ii/--"«^^^ <»«' 
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must be performed in order to obtain the observed hadronic cross section as a function of r := M'^/S. Special attention 
must be paid to the singularities in the resummed exponents ^^^' , which are situated at A = 1/2 and are related 
to the Landau pole of the perturbative coupling as- To avoid this pole as well as those in the Mellin moments of 
the PDFs related to the small- a; (Regge) singularity fa/A{x,Ho) ex a;"(l — x)^ with a < 0, we choose an integration 
contour Cn according to the principal value procedure proposed in Ref. [33| and the minimal prescription proposed 
in Ref. [3J| and define two branches 

Cat: N^C + ze^^'f' with ze[0,oo[, (99) 

where the constant C is chosen such that the singularities of the iV-moments of the PDFs lie to the left and the 
Landau pole to the right of the integration contour. While formally the angle can be chosen in the range [jt/2,tt[, 
the integral converges faster, if > tt/2. 

IV. NUMERICAL RESULTS 

We now turn to our numerical analysis of threshold resummation effects on the production of various gaugino pairs 
at the Tevatron pp-colhder {^/S — 1.96 TeV) and the LHC pp-coUider {\/S = 7—14 TeV). For the masses and widths 
of the electroweak gauge bosons, we use the current values of m^ = 91.1876 GeV and mw = 80.403 GeV. The squared 
sine of the electroweak mixing angle 

sin^^H- = 1-^ (100) 

m^ 

and the electromagnetic fine structure constant 

^/2GFml^, sin^ 9w ^Q-^^ 

can be calculated in the improved Born approximation using the world average value of Gp — 1.16637 • 10^^ GeV"^ 
for Fermi's coupling constant ■35'. The CKM-matrix is assumed to be diagonal, and the top quark mass is taken to 
be 173.1 GeV 36]. The strong coupling constant is evaluated in the one-loop and two-loop approximation for LO 
and NLO/NLL-f NLO results, respectively, with a value of A-^ corresponding to the employed LO (CTEQ6.6L1) 
and NLO (CTEQ6.6M) parton densities [37]. For the resummed and expanded contributions, the latter have been 
transformed numerically to Mellin A^-space. When we present spectra in the invariant mass M of the gaugino pair, 
we identify the unphysical scales fip — IJ-r — M with M, whereas for total cross sections we identify them with 
the average mass of the two produced gauginos. The remaining theoretical uncertainty is estimated by varying the 
common scale ^ about these central values by a factor of two up and down and the parton densities along the 22 
eigenvector directions defined by the CTEQ collaboration. 

A. Benchmark points 

The running electroweak couplings as well as the physical masses of the SUSY particles and their mixing angles are 
computed with the computer program SPheno 2.2.3 [3a], which includes a consistent calculation of the Higgs boson 
masses and all one-loop and the dominant two-loop radiative corrections in the renormalization group equations 
linking the restricted set of SUSY-breaking parameters at the gauge coupling unification scale to the complete set of 
observable SUSY masses and mixing angles at the electroweak scale. We choose the widely used minimal supergravity 
(mSUGRA) point SPSla' [20] as the benchmark for most of our numerical studies. This point has an intermediate 
value of tan/3 =10 and ^ > (favored by the rare decay 6 — > 57 and the measured anomalous magnetic moment of 
the muon), a light gaugino mass parameter of 7tt.i/2 — 250 GeV, and a slightly lower scalar mass parameter ruQ = 70 
GeV and trilinear coupling Aq = —300 GeV than the original point SPSla |^] in order to render it compatible with 
low-energy precision data, high-energy mass bounds, and the observed cold dark matter relic density. It is also similar 
to the post-WMAP point B' (mo = 60 GeV and Aq = 0) [40|, which has been adopted by the CMS collaboration as 
their first low-mass point (LMl) [4l|. In the SPSla' scenario, the Xi is the LSP with a mass of 98 GeV, the gauginos 
producing the trilepton signal have masses of Ti-i ~ m^o = 184 GeV, and the heavier gauginos, which decay mostly 
into the lighter gauginos, W and Z bosons as well as the lightest Higgs boson, have masses of m^o — 400 GeV and 
m-± ~ m^o = 415 GeV. The average squark and gluino masses are m.q ~ 550 GeV and m,g = 604 GeV. 
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TABLE II: Names, mSUGRA parameters, and physical SUSY particle masses of the benchmark points used in our numerical 
studies. 



Scenario 


mo [GeV] 


mi/2 [GeV] 


Ao [GeV] 


tan^ 


sgn(^) 


m^o [GeV] 


m-.,,„[GeV] 


m^o [GeV] 


m,,.,,o[GeV] 


m,- [GeV] 


m-g [GeV] 


SPS la' 


70 


250 


-300 


10 


+ 


98 


184 


400 


415 


550 


604 


LMl 


60 


250 





10 


+ 


96 


178 


346 


366 


550 


603 


LM7 


3000 


230 





10 


+ 


94 


176 


337 


359 


3000 


636 


LM9 


1450 


175 





50 


+ 


70 


128 


263 


284 


1480 


487 


SU2 


3550 


300 





10 


+ 


124 


229 


355 


384 


3560 


809 


SU3 


100 


300 


-300 


6 


+ 


118 


223 


465 


481 


650 


715 



Apart from the low-mass point LMl, we will also study the points LM7 and LM9, since all three points have been 
found by the CMS collaboration to lead to visible three-lepton signals. For LM7, the direct X]'^X2 production cross 
section exceeds even 70% of the total SUSY particle production cross section [4l|. The ATLAS collaboration have 
studied the direct production of gauginos at the points SU2 and SU3 with or without a jet veto (denoted JV, i.e. no 
jet in the event with transverse momentum pT > 20 GeV) in order to suppress the background from top quark pair 
production ^^. A summary of all scenarios considered here is presented in Tab. [Ill Note that none of these points 
falls into (but most of them lie relatively close to) the regions excluded by the Tevatron collaborations CDF and DO, 
which assume, however, a lower value of tan/3 = 3 and always ^o = [la]. In Ref. [14], the cross section for the 
associated production of xf and X2 ^^ been computed as a function of tan/3 and TO1/2 for mo = 200 and 1000 GeV 
and assuming Aq = and /i > 0. Unfortunately, the exact version of the renormalization group program SPheno 
used there could not be determined, and we were not able to reproduce the physical SUSY particle mass spectra of 
Ref. [13] ■ Since we also do not completely agree analytically with the coefficient function "H^j^ (M^, /x^) of Ref. [13] 
(see above), we must refrain from a direct comparison of our numerical results. 



B. Invariant mass spectra 



In Fig.|n]we present invariant mass spectra AI^dcr/dM for the production of various combinations of Xj^ and X2 with 
TO-± ~ m^o = 184 GeV in the SPSla' scenario at the Tevatron. The spectra start at M = m-± + m^o = 368 GeV 

Xi X2 -^ Xi X2 

and increase considerably from LO (blue) to NLO (grey), but much less from NLO to NLL+NLO (red). The scale 
uncertainty is considerably reduced from NLO to NLL+NLO, which indicates good convergence of the reorganized 
perturbative series. The cross section is smallest for the production of two neutralinos, since they are gaugino-like 
and couple only weakly to the s-channel Z-boson (see Tab. H]). Since the Tevatron is a pp collider, the cross sections 
are identical for X^XT ^^'^ Xi X2 (^^^ shown) pairs. The largest cross section is obtained for chargino pairs due to 
the s-channel photon contribution. Threshold resummation should be most important as M — >■ ^/s and z — >■ 1, but 
its effects on the partonic cross section are, of course, reduced in the hadronic cross section shown here by the parton 
densities, which tend to as Xa^b and z — >■ 1. Nevertheless, on close inspection one observes that the NLL-I-NLO cross 
section for two neutralinos no longer overlaps with the one at NLO for relatively large invariant masses of M ~ \/S/2. 

A similar hierarchy of the different production channels is observed in Fig. [10] for the LHC with its current center- 
of-mass energy of vS — 7 TeV and in Fig. [TT] for the LHC with its design energy of ^/S — 14 TeV. There are, however, 
two notable differences. First, the LHC is a pp collider, so that the cross section for X]^X2 exceeds the one for xrX2 
by a factor of two and becomes even larger than the one for chargino pairs. Second, the NLO band is separated by 
a wider gap from the LO band than it was the case at Tevatron, whereas the NLL-I-NLO and NLO bands overlap 
considerably more. This is, of course, due to the fact that the light gauginos are now produced further away from 
the threshold of the 7 or 14 TeV collider, so that the importance of soft-gluon resummation is reduced. However, one 
still observes a sizeable reduction of the scale uncertainty from NLO to NLL-I-NLO. 

Heavier gaugino pairs can only be produced with sizeable cross sections at the LHC. We therefore show in Figs. 
1121 andlTSlthe invariant mass spectra M^da/dAI for the production of various combinations of X3 4 and X2 at the 
LHC with ^/S ^7 TeV and 14 TeV and with m^o = 400 GeV and m~± ~ m^o = 415 GeV in the SPSla' scenario. 

X3 A2 '^4 

The spectra start at M ~ 800 — 830 GeV, and their magnitudes are considerably smaller than in the light gaugino 
case. However, they are now of comparable size for neutralino and chargino pairs due to the fact that the dominantly 
higgsino x-j and Xa now have sizeable couplings to the s-channel Z-boson (see Tab.|I]). The associated production of a 
neutralino and a chargino is again much larger for the positive chargino eigenstate than for its negative counterpart. 
The cross sections for X3X2 pairs are very similar to those for X4X2 pairs and therefore not shown. Higgsino-like 
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FIG. 9: Invariant mass spectra for the production of various light gaugino pairs at the Tevatron in the SPSla' scenario and in 
the LO (blue), NLO (grey) and NLL+NLO (red) approximation. The corresponding scale uncertainties are represented by the 
band widths. 



neutralinos and charginos with large s-channel contributions arc produced as iS-waves, so that the invariant mass 
spectra rise more steeply at low M than P-wave produced gaugino-like neutralinos and charginos. 

From Figs. l9lfT3l the impact of threshold resummation effects is difficult to estimate. We therefore present in Fig. 
[14] the relative size 



^NLL ^ 



^^NLL+NLO 
dcrNLO 



(102) 



of the NLL+NLO prediction with respect to the NLO prediction. As one expects, the correction is larger at the 
Tevatron with its lower center-of-mass energy (top left) than at the LHC (top right and bottom) and increases with 
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FIG. 10: Same as Fig. [9] for the LHC with its current center-of-mass energy of yS = 7 TeV. 



the invariant mass. The relatively small differences among the i^^'^'^-factors for neutralino pair production and the 
channels involving at least one chargino can be traced to the fact that the former receives most of its contributions from 
t- and M-channel squark exchanges, which are more sensible to strong corrections than the exchanges of electroweak 
bosons in the s-channel. 

The JC^'^'^-factors for the production of heavy gaugino pairs at the LHC with ^/S — 7 TeV (14 TeV) are presented 
in Fig. llSl left (right). They are larger than their counterparts for light gauginos in Fig.ll4| since the gaugino masses as 
well as the invariant masses M are now closer to the hadronic center-of-mass energies. In addition, the result for the 
X3X4 channel differs no longer substantially from the other channels, since the heavy neutralinos are now higgsino-like 
and their production is now also dominated by the s-channel exchange of a weak gauge boson. 
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FIG. 11: Same as Fig. [9] for the LHC with its design center-of-mass energy of v5 = 14 TeV. 



C. Total cross sections 

The stability of the perturbative series and its reorganization is traditionally checked by varying the factorization 
and renormahzation scales ^f and f^u about a central value hq. We therefore present now the total cross sections for 
the production of light gaugino pairs at the Tevatron f Fig. [TBI) and at the LHC with v^ = 7 TeV fFig. [T7|) as a function 
of the ratio hf,r/ij-o, where the central scale /iq is now chosen to be the average mass of the produced gaugino pair. 
The LO prediction (blue, dot-dashed) of the electroweak processes under consideration is, of course, independent 
of the renormalization scale ^r (right part of the figures), whereas the NLO prediction (black, dashed) depends 
inversely on the logarithm of fj,R through the strong coupling asipn). At NLL accuracy (red, full), the resummed 
soft corrections attenuate this dependence and introduce a plateau region, so that the prediction is stabilized. The 
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FIG. 12: Same as Fig.|9]for the production of heavy gaugino pairs at the LHC with its current center-of-mass energy of VS = 7 
TeV. 



factorization scale /i^ (central part of the figures) enters the hadronic cross section already at LO through the largely 
logarithmic dependence of the PDFs, which is then attenuated by the factorization of initial-state singularities at 
NLO and further at NLL accuracy. In all cases, the resulting total NLL+NLO prediction is thus much less dependent 
on the common scale /zf = Mfi = M (l^ft part of the figures) than the LO and NLO estimates. 

In Tab. IIIII we present the total cross sections for the trilepton channel in the SPSla' scenario at the Tevatron 
{yfS = 1.96 TeV) and LHC {\fS = 7, 10 and 14 TeV). Besides the central values (in fb) at LO, NLO and NLL+NLO, 
we also present the scale and PDF uncertainties. The former are estimated as described above by a common variation 
of the renormalization and factorization scales by a factor of two about the average mass of the two gauginos, the 



21 



> 

o 

s 



b 



> 

a 

s 



b 



PP ^ X3X4 



-^-^ atVS = U TeV 



35 
30 
25 
20 
15 
10 
5 






iNLL+NLO 1 1 


/ 


^ NLO r;::;] - 


/ 


\. LO r"'\ - 


//- 




/ '' ' 


^^ N vV 


"1* 


"^^^\^ \v 


If 


^^^^\\v 


f 


^ x^ ^ ^s. 




''•■^^■■-. ^^ 








^■^■-'^^S^ 


^ 


1 1 1 1 1 



800 1000 1200 1400 1600 1800 2000 

M [ GeV ] 

PP -^ xtxl at v^ = 14 TeV 




800 1000 1200 1400 1600 1800 2000 

M [ GeV ] 



> 

o 

s 



b 



PP -> X4X2" at v^ = 14 TeV 








800 1000 1200 1400 1600 1800 2000 

M [ GeV ] 

PP -^ xtX2 at v^ = 14 TeV 




800 1000 1200 1400 1600 1800 2000 

M [ GeV ] 



FIG. 13: Same as Fig.[9]for the production of heavy gaugino pairs at the LHC with its design center-of-mass energy of v5 = 14 
TeV. 
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along the 22 eigenvector directions defined by the CTEQ collaboration. Since these are available only for the NLO 
fit CTEQ6.6M, but not for the LO fit CTEQ6.6L1, we do not present a PDF uncertainty for the LO prediction. 
Furthermore, the same PDF set enters at NLO and NLL+NLO, so that the PDF uncertainties for these two predictions 
coincide. The most important result is again the considerable reduction of the scale uncertainty from LO to NLO 
and then to NLL+NLO. The total cross sections increase with the available center-of-mass energy due to the higher 
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which agrees with the cross sections given in Tab. |TTT] surprisingly well. 

In Tab.|IV]we fix the LHC center-of-mass energy to its design value of v^ — 14 TeV and show the total production 
cross sections for light and heavy gaugino pairs in LO, NLO and NLL+NLO together with the corresponding theoretical 
uncertainties. As it was already mentioned above, the cross section for the higgsino-like XsXl pairs is about as large as 



23 



pp 



{VS = 7 TeV) - SPSla' 



pp 



{^/S = 14 TeV) - SPSla' 






1.065 

1.06 
1.055 

1.05 
1.045 

1.04 
1.035 

1.03 
1.025 



^0 ^0 






xixi -- 


.... 


/ 
/ 


xlx2 






- __|_ _o 




/ , 


X2X4 - 






. 1- ~ 


/ 


y J 


" X2 X2 - 




y^ 








- 


/'' •^' y^ 


..■■■■■ 




f* x/ 




y' 






y'--^'^ 




" 












































jT y 






1 


1 1 


1 



ft 




900 100011001200130014001500 

M [ GeV ] 



900 100011001200130014001500 

M [ GeV ] 



FIG. 15: Same as Fig. [H for heavy gaugino pairs at the LHC with sfS = 1 TeV (left) and ^ = 14 TeV (right). 

TABLE III: Total cross sections for the production of Xi^X2 pairs in the SPSla' scenario at different hadron colliders and 
center-of-mass energies in the LO, NLO and NLL+NLO approximation, together with the corresponding scale and PDF 
uncertainties. 
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the one for X2 X2^ pairs, and it is in fact not much smaher than the one for the considerably lighter gaugino-like X2X2 
pairs. In general, the heavy gaugino cross sections are, however, significantly smaller than those for light gauginos. 

In Tab.FVl we present finally total cross sections for the trilepton channel in our different benchmark scenarios at the 
LHC with its current center-of-mass energy of VS = 7 TeV. Since the masses of Xi and X2 &re always rather similar, 
one expects also similar total cross sections. This is indeed confirmed by the VS — 7 TeV results in Tab. [TTTland the 
numbers in Tab. |V]with the notable exceptions of LM7 and LM9, where the cross section is about a factor of two and 
one order of magnitude larger than for the other points, respectively. This is partly due to the lower gaugino masses 
at LM9 and partly to the much heavier squark masses, which suppress the t- and w-channels and thus their destructive 
interference with the s-channel amplitudes. The additional jet veto (JV), i.e. the rejection of events containing jets 
with transverse momentum pT > 20 GeV, envisaged by the ATLAS collaboration to suppress the background from 
top quark pair production, has obviously no consequences at LO, since gauginos are exclusively produced at this order. 
An additional quark or gluon can only be present at NLO or NLL-fNLO, and restricting its transverse momentum to 
low values reduces the total cross sections slightly with respect to the unrestricted predictions. The small reduction 
of the signal cross section in combination with a large reduction of the background should therefore indeed lead to a 
much better significance. 



CONCLUSIONS 



In summary, we presented in this paper a complete analysis of threshold resummation effects on direct light and 
heavy gaugino pair production at the Tevatron and at the LHC with its current, intermediate, and design center- 
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FIG. 16: Total cross sections for the production of neutralino (top left), chargino-neutralino (top right) and chargino pairs 
(bottom) at the Tevatron with ^fS = 1.96 TeV in the LO (blue, dot-dashed), NLO (black, dashed) and NLL+NLO (red, full) 
approximation. 



of-mass energies. We confirmed the existing NLO calculation and extended it to include also mixing effects for 
intermediate squarks. Soft gluon radiation in the threshold region was resummed at leading and next-to-leading 
logarithmic accuracy into a Sudakov exponential, and the full SUSY-QCD corrections were retained in the finite 
coefficient function. This allowed us to correctly match the resummed cross section at NLL accuracy to the NLO 
perturbative cross section. Universal subleading logarithms coming from the splitting of initial partons were resummed 
in full matrix form, i.e. including also non-diagonal splittings. We found that resummation increased the invariant 
mass spectra and total cross sections only slightly, but stabilized the reorganized perturbative series considerably with 
respect to the fixed-order calculation. For future reference, we presented total cross sections with the corresponding 
theoretical errors from scale and PDF variations at LO, NLO and NLL-I-NLO in tabular form for several commonly 
used SUSY benchmark points, light and heavy gaugino pairs, and various hadron collider types and energies. 
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FIG. 17: Same as Fig. [16] for the LHC with its current center-of-mass energy of VS — 7 TeV. 



Appendix A: Quark and squark couplings to weak gauge bosons and gauginos 

The coupling strengths of the electroweak gauge bosons Z and W to left- and right-handed (s)quarks q (q) with 
weak isospin T^ and fractional electric charge Cq are given by 



{LqqZ,Rqqz} 
{Lq^qjZ,Rq,qjZ} 



-^n-eqxw). 



2cw 



{Ludw,Rudw} — { — 
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2cw 



iT^-eqXw){RiR%RhRT2} ^ 
.0}, 



{^UidjW^^u.djw} - {^^^-Ril^jl' 0}' 



(Al) 



sin 6w is the squared sine of the electroweak mixing angle, R!^: are the elements of the 
rotation matrices diagonalizing the up- and down-type squark mass matrices, and Vud are the elements of the CKM 



where xw = 1 — c- 



w — ''W 
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TABLE IV: Total cross sections for the production of various gaugino pairs in the SPSla' scenario at the LHC with its design 
center-of-mass energy of \/S = 14 TeV. The central predictions are given at LO, NLO and NLL+NLO together with the 
corresponding scale and PDF uncertainties. 
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TABLE V: Total cross sections for the production of X+X2 pairs at the LHC with its current center-of-mass energy of ^fS — 7 
TeV for different SUSY benchmark points. The central predictions are given at LO, NLO and NLL+NLO together with the 
corresponding scale and PDF uncertainties. 
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matrix. Their SUSY counterparts, i.e. the squark-quark-gaugino couphngs, are given by 
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where the matrices iV, t/ and V diagonahze the neutral and charged gaugino-higgsino mass matrices, m^d are the up- 
and down-type quark masses, and mw is the mass of the VF-boson. All other couplings vanish due to electromagnetic 
charge conservation [Ij]. 
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